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Abstract 

We introduce the standard distribution width of fitness to characterize the 
global and individual features of a ecosystem in the Bak-Sneppen evolution 
model. Through tracking this quantity in evolution, a different hierarchy of 
avalanche dynamics, wq avalanche is observed. The corresponding gap equa- 
tion and the self-organized threshold w c are obtained. The critical exponents 
r, 7 and p, which describe the behavior of the avalanche size distribution, 
the average avalanche size and the relaxation to attractor, respectively, are 
calculated with numerical simulation. The exact master equation and 7 equa- 
tion are derived. And the scaling relations are established among the critical 

exponents of this new avalanche. 
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I. INTRODUCTION 



The biological evolution theory of Darwin is well known to all. The biological evolution 
can be divided into two different types, micro-evolution and macro-evolution, which focus on 
the dynamics of evolution at the level of a single population and of all species in ecosystem 
respectively. Up to now, many scaling laws have been found in macroevolutionary data. The 
distribution N(m) of extinction sizes for families m decreases with m according to a power- 
law: N{m) oc m~ a 0. The distribution N(t) of genera lifetimes t follows N(t) oc t~ k 0. 
The distribution of the number of genera N(g) with g species can be fitted with N(g) oc g~^ 
||. The disappearance of species or genera is similar to radioactive decay, in which the 
amount of the "original element" decays exponentially with time. This phenomenon has 
been named as Red Queen effect [|J. The extinction events exhibit punctuated patterns 0. 
The punctuated equilibrium can be interpreted as a metastable state of the dynamics of that 
single species ||. If an ecosystem is in equilibrium, every species is stable. However, if one 
of them fluctuates with time, it will lead to the destabilization of its neighboring species. 
Thus, the avalanche fluctuation will sweep through the evolution. 

The existence of scaling laws in macro-evolution indicates that the dynamics of large-scale 
evolution is the result of a self-organized critical process |7]^9|]. Based upon this view, Bak 
and Sneppen introduced a toy model for species evolution pTOfl . It is defined and simulated 
as follows: (i) L species are arranged on a one-dimensional lattice with periodic boundary 
conditions, (ii) Random numbers fi, which are chosen from a uniform distribution between 
and 1, P(f), are independently assigned to each site (species), (iii) At each time step, the 
extremal site, i.e., the one with the smallest random number / m i n , and its nearest neighbors 
are mutated by assigning new random numbers in P(f). As the system evolves it turns 
out that the species with fitness significantly above the self-organized threshold f c (equal to 



0.66702 [|Tl|]) will never mutate, unless its neighbors mutate. After enough mutations have 
occurred, the ecosystem reaches a statistically stationary state in which the distribution for 
the fitness of all species is statistically stationary. 
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Since the Bak-Sneppen model was denned, its punctuated patterns and avalanche dynam- 
ics have been investigated extensively. Through observing the fluctuations of the smallest 
fitness, the model shows punctuated equilibrium behavior, i.e., it self-organizes into a critical 
state with intermittent coevolutionary avalanches of all size ||10|| . Paczuski, Maslov and Bak 
studied the f avalanche (PMB avalanche), forward avalanche and backward avalanche and 
derived the exact scaling equations, they also established the relations among the critical ex- 
ponents [fOll . Recently, by monitoring the variations of the average fitness, the / avalanche 



(LC avalanche), is found, and the exact equations and scaling relations are obtained fl2|-|l4 
As we know, the avalanche dynamics is a kind of macroscopic phenomenon in driven dis- 
sipative systems. The detailed dynamics of these systems sensitively depend on the initial 
configurations. But the distribution of avalanches, i.e., scaling law, does not depend on 
such details due to the universality of complexity (T^|. In the existing works, analyzing the 



smallest fitness only directly connects with the feature of individuals, but it does not directly 
represent the global feature; discussing the average fitness only directly connects with the 
global dynamics, but it does not directly manifest the difference in individuals. Can we find 
a new quantity, which can describe the global feature and the individual difference simulta- 
neously? In this paper, we define such a quantity and observe the corresponding avalanche 
dynamics with this quantity. 

In the next section, the standard distribution width of fitness is defined, and its gap 
equation is obtained. In Sec. Ill, the dynamics of avalanche is analyzed in detail, the 
exact mast equation is derived, and the scaling relations are established among the critical 
exponents. The last section is a brief summary. 

II. STANDARD DISTRIBUTION WIDTH OF FITNESS AND ITS GAP 

EQUATION 

In the Bak-Sneppen model, the fitness of the i-th species is denote as fi (i — 1, 2, • • 
■,L), with L being the total number of species. The fitness represents the population or 
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the living capability of the species, large fitness means immense population or great living 
capability, and vice versa. To describe the average population or average living capability 
of an ecosystem, Li and Cai defined the average fitness as [p^ -Tl4]| 

f av {s) = \ir,m. (i) 

i=l 

where fi(s) is the fitness of the i-th species at time s. Based upon the average fitness, we 
define the standard distribution width of fitness as 



w(s) 



i 



^E(/.w-/»w) 2 . (2) 

i=i 

From the above definition, we know that the standard distribution width of fitness simulta- 
neously manifests both global dynamics and individual difference. At the beginning of the 
evolution, i.e., s = 0, the fitnesses are uniformly distributed in (0, 1). Thus for an infinite- 
size ecosystem (L — > oo), the initial average fitness / a „(0) equals 0.5, the initial standard 
distribution width w(0) is \/3/6. But for a finite-size ecosystem, f av (0) and w(0) fluctuate 
around the above values due to the finite size. Because of the absence of mutations or up- 
dates, w(0) does not show the correlation between species, only when the evolution goes on, 
such correlation is contained in w(s) gradually. 

As time increases, f a v{s) shows an increasing tendency ||12|| , whereas, w(s) shows a de- 



creasing tendency (see the bottom of Fig. 1). After a long transience, i.e., s — > oo, the ecosys- 
tem self-organizes to a critical state, and, f av (s) and w(s) approach to different thresholds. 
When s — > oo, L — > oo, almost all species have fitnesses above the threshold f c = 0.66702 
IT!] , and the fitnesses are uniformly distributed in (f c , 1). With the definitions of f av {s) and 



w(s), we obtain the relation among the thresholds 

^=^(Wc) = ^(l-/a,c). (3) 

Here, w c and f avc are thresholds of w(s) and f av {s) respectively. Thus, one gets w c = 0.09623 
and favc = 0.83351. For the finite-size ecosystem, the self-organized critical thresholds will 
fluctuate around the above values, and this fluctuation decreases with the increase of L. 



Similar to the definition in Refs. fi~l]l2 |, we define the gap of w(s) as its envelope: 



, w(0) , 8 = 0, 

ir(.s) = < (4) 

mm(w(s)\s' E [0, s]) , s > 0. 

The above definition means that the current gap of w(s) is the most minimum of all w(s A ) 
(0 < s' < s). In the top of Fig. 1, we show the gap W(s) as a function of s. With the 
increase of time, W(s) gradually approaches close to the threshold. By definition, when the 
gap falls to its next lower value, the separate instances are separated by avalanches. The 
avalanches correspond to plateaus in W(s) during which w(s) > W(s). A new avalanche 
is initiated when the gap falls, and the old avalanche terminates at the same time. Once a 
avalanche is over it will never affect the behavior of any subsequent avalanche. When the 
gap decreases, the probability for w(s) to jump above W(s) increases, thus longer and longer 
avalanches happen typically. Finally, the ecosystem falls into statistically stationary state. 

Following Refs. we derive the gap equation of W(s). If the gap is changed from 

W(s) to W(s) + AW, the average number of avalanches occurred is N av = — AWL/(W(s) — 
w c ). Thus, when L — > oo, the average number of time steps required for increasing the 
gap from W(s) to W{s) + AW is As = (S) w(s) N av = - (S) w(s) AWL/(W(s) - w c ), where 
(S) w(s) ls ^ ne average size of avalanche of the current plateau of the gap. So the differential 
equation for W(s) is 

dW(s) , AW Wis) -w c 

— = hm — — = — — — . (5) 

ds AWAs^o As (S) W ,,L 

Eq.(5) describes the relaxation to attractor. If the law of {S) w ^ is obtained, one can derive 
the law of relaxation to attractor from the above equation. 

III. AVALANCHE DYNAMICS 

We know that all self-organized critical systems will exhibit a power-law avalanche 
dynamics. This has been confirmed by a lot of models with self-organized criticality 
|ll3|- |i^[i6H 22|] . In this section, based on the definition of PMB avalanche IJ] and LC 



avalanche fL2| , the definition of a new avalanche, Wo avalanche, is presented. Then its 
critical exponents are calculated with numerical simulation and its exact master equation is 
derived. Lastly, the scaling relations are established. 



The definition for wq avalanche is as follows. Similar to those used in Refs. fTTyizfl , we 
introduce an auxiliary parameter w , where w(0) > w > w c . Suppose that at time step 
s\, the current standard distribution width w(si) is less than wq. If the next standard 
distribution width w(si + 1) is larger than wq, a creation-annihilation branching process is 
initiated. The avalanche will continue to run unless w(s) becomes less than wq. This means 
that, at time s, if all w(s) > wq for (s\ < s' < s — 1), the current avalanche continues. 
In terms of the above definition, the size of an avalanche is the number of time steps of 
subsequent punctuations above w . If the first appearance of w(s) < w Q after time si occurs 
at time S\ + S, the size of the current avalanche is S. 

Clearly, the above definition ensures the hierarchical structure, i.e., larger avalanches 
consist of smaller ones. With the decreasing of Wq, smaller avalanches combine into larger 
ones; if Wq is set as w c , the infinite-size avalanche appears. In the other case, with the 
increasing of w$, larger avalanches split into smaller ones. Thus, the cutoff effects is un- 
avoidable when Wq is not chosen as w c . Nevertheless, the same scaling laws can be obtained 
when wq is very close to the threshold w c . For the Bak-Sneppen model with 200 species, 
choosing Wq as 0.1100, we find by numerical simulation that the distribution of avalanche 
size follows the power law 

P{S) oc S~ T , (6) 

where P(S) is the probability of avalanche of size S, r equals 1.63 ± 0.06 (see Fig. 2). 
Because of the cutoff effects, the average avalanche size varies when wq changes. Similar to 
those in Refs. [Tl]|T^,[14]] , our numerical results show that the average size of Wq avalanche 



obeys 

(S) wo cx(w -w c )~^ (7) 
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Here, w c is the threshold, the exponent 7 equals 2.57±0.05 (see Fig. 3). This law means that 
the divergence of the average avalanche size satisfies power law. Our numerical simulation 
also shows that the relaxation to attractor abides 

(W - w c ) a s' p . (8) 

Here, s is the time steps, the exponent p equals 0.39 ± 0.03. Obviously, the value of p is 
very close to I/7. 

So far, we have found some scaling laws of the wq avalanche by numerically simulations. 
Now, we will give some exact results. To describe the cascade process of smaller avalanches 
combining into larger ones when the parameter wq is changed, an exact master equation 
is derived as follows. Denoting the probability of w avalanche of size S as P(S,w ), when 
parameter changes from u> into w + Aw , some w avalanches merge to larger w + Aw 
avalanches, and some small avalanches merge into Wq avalanche. Thus, at the same time, 
the probability of a avalanche with a given size flows in and out. Because the termination 
of wo avalanche is uncorrected, the probability of a wo avalanche merging to wo + Au> 
avalanche is proportional to —Awq/{wq — w c ) [|1^,|13|]. Thus the probability which flows out 
from P(S, Wq) can be written as 

A out P = -XP(S, w Q )Aw /(w - O, (9) 

where A is a constant. The probability which flows in P(S,wq) is 

5-1 

A m P = -\J2 P ^ w ^ Pi - S ~ 5l ' ™o) ^ /(w - w c ). (10) 

Sx=l 

When wo approaches w c and Au> approaches zero, we obtain the following master equation 

S—l 

_ ( Wo _ Wc ) dP £> w °) = _ X P(S, w ) + Xj2 p ( S ^ w °) p ( S - s u wo). (11) 

W ° Si=l 

As pointed out in the previous section, with the increasing of time, the standard distribution 
width shows a decreasing tendency. So, when wq decreases, the first term on the right side 
of the master equation reflects the loss of S'-size avalanche, and the second term depicts the 
gain of S'-size avalanche. 
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Similar to Ref. |n|, after defining a new variable u = — \n(wo — w c ), the master equation 
becomes 

s-i 



9P i S,U) = ~XP(S, u) + \J2 P(S 1 ,u)P(S - St, u). (12) 
Through performing a Laplacian transformation for P(S,u), i.e., 

oo 

p(a,u) = J2P(S,u)e- aS , (13) 
5=1 

thus the master equation reads as 

^^^Xp(a,uMa,u)-l]. (14) 

oo 

When a = 0, p(a, u) = P(S,u)e~ aS = 1, this is the normalization of P(S,u). While 

s=i 

a > 0, one gets < p(a, u) < 1. Expanding both sides of Eq.(14) as Taylor series throughout 
the neighborhood of the point a = 0, one can obtain 

^(S) U -^(S\ + ^(S\...}, (15) 

oo 

where \S k ) = S k P(S, u). Comparing the coefficients of different powers of a in Eq.(15), 

5=1 

one gains an infinite series of exact equations. The first exact equations can be expressed as 

i^S) u = X. (16) 
Replacing u as — ln(w7 — w c ), we obtain the 7 equation 

ln(5L^--A_. (17) 



dw w ° w - w c 

Noting the scaling law (S) w oc (w — w c )~ 7 , inserting it into the above equation, one can 
easily gain 

7 = A. (18) 
8 



By far, we have obtained three critical exponents of the wq avalanche: r, 7 and p, which 
describe the behavior of the avalanche size distribution, the average avalanche size, and the 



relaxation to attractor, respectively. As listed in Ref. |TT[, there are three other critical 
exponents: D, a and v. The avalanche dimension D is defined by the following scaling 
relation 

S oc R D . (19) 

where S is the avalanche size (temporal duration), R is the spatial extent. The a and v are 
defined by the following scaling laws of cutoff 

S co cx (w -w c )- 1/tT , (20) 

R co cx (w - w c )~ v . (21) 

The symbols S co and R co are the cutoff of the avalanche size and that of spatial extent of 
the avalanche. 

Although there are six critical exponents, only part of them are independent. After 
choosing r and 7 as the independent exponents, the other ones can be expressed as functions 
of them. Inserting the scaling law (S) w oc (wo — w c )~ 7 into the gap equation (Eq.(5)), and 
noting the scaling law (wq — w c ) oc s~ p , one can find 

p = 7~ 1 . (22) 



Integrating the equation {S) w = $ l c ° SP(S,w )dS, it is easy to obtain 

/Sco rSco 
SP(S,w )dS = J S 1 - T dS<xS 2 - T <x(w -w c )-^. (23) 

Comparing with (S) w oc (wq — w c )~' y , one can get 

Noting that the average number of the covered sites (n cov ) Wo covered by a w avalanche 
scales near the critical point as JXTJ] (n cov ) w cc R cov cc R oc (w — tt; c ) _1 and S cc R D , so 



' w 



) W0 a 5l/D ' integrating the equation (n cov (S)) wo = J x co n cov (S)P(S,w )dS gives 
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(S)) Wo cx S l / D - T dS cx (w - » c )-^. (25) 



Thus, the avalanche dimension D can be written as 

D = = J . (26) 

r + a-1 7(r-l) + 2-r 

Noting S oc R D and i? co a (w — w c ) _tJ , one can immediately obtain 

S C0 <xR° <x(w Q -w c y vD . (27) 
Comparing with S co oc (u> — w c ) _1 / CT , it is easy to know 

Up to now, the scaling relations among the six critical exponents have been successfully 
established. The existence of two independent critical exponents means that there are two 
kinds uncorrelated critical behavior in wq avalanche. 

IV. SUMMARY 

In conclusion, by defining a new quantity, the standard distribution width of fitness, a 
new avalanche, wo avalanche, is observed in the evolution of Bak-Sneppen model. We declare 
this avalanche as a different hierarchy of avalanche, because firstly, the observed quantity is 



different from those of Refs. |10|,|TT|| and Refs. [12-14]; and secondly, the values of the critical 



exponents are also different from those found in Refs. flQjTT[| and Refs. |12|-^4[. 

From the definition of the standard distribution width of fitness, we easily obtain the 
corresponding gap equation and the self-organized threshold. Apparently, with this defini- 
tion, the global feature and the difference among individuals are simultaneously described 
by the standard distribution width. So, for the evolutionary distribution with changeless 
or small changed average values, the standard distribution width is a good quantity to ob- 
serve the evolutionary dynamics. According to the definition of PMB avalanche and LC 
avalanche, the definition of the wo avalanche is presented. With numerical simulation, the 
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critical exponents r, 7 and p are obtained. Then the exact master equation is derived. From 
the master equation, the 7 equation is immediately obtained. Finally, combining all scaling 
laws, the scaling relations are successfully established among the critical exponents. With 
these relations, we find there are only two uncorrected critical behavior in the wo avalanche, 
for there exist only two independent critical exponents. 
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FIGURES 

FIG. 1. The standard distribution width and its gap. The above row corresponds to the 
variation of the gap W(s), the bottom row corresponds to the fluctuation of w(s). 



FIG. 2. The power-law of the avalanche size distribution. 



FIG. 3. The power-law of the divergence of the average avalanche size. 
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